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One  of  the  most  important  considerations  in  the  design  of  components  of  gas  turbine  engines  is  the 
reduction  of  the  amplitudes  of  steady-state  vibrations  of  some  components  of  the  engine  during  operation. 
Large  repetitions  of  high-amplitude  vibrations  can  result  in  catastrophic  fatigue  failures  of  such  components. 
This  problem  has  become  more  critical  with  the  current  tendency  to  design  engines  to  avoid  material  inelastic 
behavior,  because  this  practice  produces  also  a  significant  reduction  of  the  amount  of  energy  dissipated  by 
means  of  viscous  or  hysteretic  damping.  Therefore,  to  keep  the  vibration  amplitudes  within  acceptable  limits, 
the  deeigners  of  gas  turbine  engines  are  using  alternative  sources  of  energy  dissipation.  The  most  widely 
accepted  solution  is  the  inclusion  of  especially  designed  friction  interfaces,  which  have  shown  good 
performance  in  laboratory  tests  as  well  as  in  airplane  operation. 

The  study  of  steady-state  vibrations  of  structures  with  so-called  frictional  dampers  has  attracted  the 
attention  of  many  researchers.  The  first  analytical  solution  for  a  single-degree-of-freedom  system  with  a 
Coulomb  frictional  joint  was  presented  by  Den  Hartog  in  1931  [1,2].  The  most  relevant  contributions  until  1979 
are  discussed  by  Plunkett  [3].  More  recently,  further  analytical  and  experimental  research  has  been  reported 
by,  among  others  Griffin  and  coworkers,  [4,5,6]  Muzynska  et  al,  [7]  Srinivasan  et  al,  [8]  Dominic  et  al,  [9]  and 
Soni  and  Bogner  [10]. 

Most  analytical  studies  are  limited  to  cases  with  only  one  friction  joint,  mainly  due  to  the  difficulties  of 
solving  the  nonlinear  differential  equations  which  result  from  the  mathematical  models.  The  available 
mathematical  models  and  methods  of  solution  cannot  be  easily  extended  to  problems  with  several  friction 
joints,  and  their  results  do  not  explain  completely  all  the  experimentally  observed  phenomena.  For  these 
reasons,  a  more  efficient  and  rational  utilization  of  friction  dampers  still  requires  additional  experimental  and 
analytical  research. 

This  report  is  the  result  of  initial  work  aiming  at  the  development  of  a  simple,  approximate,  yet  accurate, 
methodology  to  study  the  steady-state  response  of  multi-degree-of-freedom  structures  containing  friction 
interfaces,  but  which  are  otherwise  linear.  As  a  first  step  in  pursuing  this  goal  we  consider  the  steady-state 
vibrations  of  the  two-body  system  shown  in  Fig.  1,  which  is  the  simplest  case  of  linear  substructures 
connected  by  friction  joints.  The  springs  and  the  dampers  are  considered  to  be  linear,  and  the  friction  joint 
connecting  the  bodies  is  of  Coulomb  type. 

In  studying  this  problem  our  main  objectives  are:  a)  to  find  an  exact  analytical  solution;  b)  to  assess  the 
accuracy  and  applicability  of  some  approximate  methods  of  solution;  and  c)  to  understand  the  role  of  the 
friction  joint  in  the  dynamic  behavior  of  the  system. 


The  first  part  of  the  report  deals  with  the  formulation  of  the  differential  equations  governing  the  motion  of 
the  system.  Then,  an  exact  analytical  method  of  solution  for  these  equations  is  presented.  Exact  results  for 
selected  systems  are  included,  to  understand  the  role  of  different  parameters  in  the  steady-state  response. 

In  a  subsequent  section  a  complex  variable  version  of  the  harmonic  balance  method  to  obtain  approximate 
solutions  is  presented.  The  advantages  and  limitations  of  this  solution  method  are  discussed.  In  particular, 
some  modifications  to  improve  the  accuracy  of  the  method  when  the  joint  is  only  partially  sliding  are 
proposed. 

Selected  numerical  results  are  obtained  to  compare  exact  and  approximate  solutions.  The  possibility  of 
extending  the  different  mathematical  models  and  methods  to  more  general  problems  is  also  examined. 

2  EQUATIONS  OF  MOTION 

The  system  under  study  is  depicted  in  Fig.  1.  M1  and  M2  are  the  masses  of  the  bodies;  K1  and  K2  are  the 
stiffnesses  of  the  linear  springs;  C,  and  C2  are  the  constants  of  the  viscous  (or  if  it  is  the  case,  structural) 
dampers;  p  and  N  are  the  friction  coefficient  and  the  constant  normal  force  in  the  Coulomb  friction  joint.  The 
forces  applied  on  M,  and  M2  are  harmonic  with  the  same  frequency  W  and  with  amplitudes  Q,  and  Q2, 
respectively;  p  is  a  prescribed  phase  angle  of  the  second  force  with  respect  to  the  first  one. 

We  will  first  formulate  the  equations  of  motion  of  the  system  for  the  two  possible  states  of  the  friction  joint 
(slipping  and  stuck),  and  the  mathematical  representation  of  the  conditions  producing  changes  from  one  to 
the  other. 

Let  X,  and  X2  be  the  displacements  of  the  corresponding  masses,  and  F  the  tangential  force  in  the  joint. 
Then,  Newton's  second  law  yields  the  following  equations  of  motion  for  the  masses: 

M,X,  +  C,X,  +K,X,  =Q1cosWT-F  (1) 

M2X2  +  C2X2  +  K2X2  =  Q2cos(WT  +  P)  +  F  (2) 

where  dots  denote  derivatives  with  respect  to  time. 

When  the  joint  is  sliding,  F  must  satisfy  the  following  conditions: 

a)  |F|  =  pN 

b)  Its  sign  is  opposite  to  that  of  the  relative  velocity,  XrX2. 

These  conditions  reflect  the  nonlinear  nature  of  the  problem,  because  both  X,  and  X2  depend  on  F,  which,  in 
turn,  depends  on  their  derivatives. 
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In  order  to  identify  the  parameters  governing  the  dynamic  response  of  the  system  and  to  obtain  results 
applicable  to  a  wider  range  of  cases,  it  is  convenient  to  express  the  above  equations  in  a  dimensionless  form. 
For  this  purpose,  define  the  following  constants: 
m  =  +  M2) 

k  =  Kg/tK,  +  K2) 

T0  =  f(M,  +  M2)/(K,  +  K2)]1/2 
w  =  WT0 

Xs  =  Q1/(K,  +  K2) 
f  ■•=  F/Q, 

<0  =  pN/Q, 

C,  =  Cj/[2  (Kj  Mj)1/2]  ,  j  =  1 ,2 
=(K)/Mj),/2  .j=1,2 

u  =  Q^Q,  (3) 

and  the  following  dimensionless  variables: 

Xj-X/X,  .i  -  1.2 

t  =  T/To  (4) 

With  the  above  definitions  simple  algebraic  manipulations  transform  (1)  and  (2)  into 
(1-m)  x,  +  2C1lF21x1  +  (1-k)  x,  =  -f0  +  coscot  (5) 

m  x2  +  2^2X2  +  k  x2  =  f0  +  acos((ot  +  P)  (6) 

These  are  the  nondimensional  equations  for  the  system  of  Fig.  1 .  The  conditions  to  be  fulfilled  by  the  force  f 
when  the  joint  is  sliding  become 
!f|  =  <o 

sign(f)  =  -sign(x1  -  x2) 

When  the  joint  is  stuck  the  difference  between  the  displacements  of  the  masses  is  equal  to  some  constant 
value  d,  i.  e. 

x,=x2  +  d  (7) 

As  a  consequence,  the  time  derivatives  of  x,  and  x2  are  the  same  and  equations  (5)  and  (6)  reduce  to 

x,  +  2(C1Q1  +  +  x1  =  cos(o)t)  +  acos(cot+P)  +  kd  (8) 

This  equation  holds  while  the  amplitude  of  the  tangential  force  in  the  joint  remains  below  f0. 
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Exact  analytical  solutions  for  the  differential  equations  of  motion  are  desirable  because  they  describe  all  the 
characteristics  of  the  system  response.  Hence,  the  significance  of  all  the  physical  parameters  can  be 
precisely  determined.  This  eventually  will  provide  a  basis  for  proposing  approximate  methods  of  solution  and 
also  the  means  to  assess  the  accuracy  and  limitations  of  those  methods. 

In  our  case,  the  difficulties  to  obtain  analytical  solutions  arise  from  the  fact  that  it  is  neccesary  to  know  at 
any  time  whether  the  friction  joint  is  stuck  or  sliding  in  order  to  use  the  corresponding  equation  of  motion. 
However,  this  information  is  part  of  the  solution  rather  than  part  of  the  data.  In  generd,  an  unknown  number 
of  incursions  in  both  conditions  can  take  place  in  a  semiperiod  of  vibration. 

In  the  analytical  method  presented  in  this  section  it  is  assumed  that  the  joint  is  either  always  sliding  or  it 
only  experiences  one  slip-to-stuck  change  during  a  semiperiod.  The  first  step  of  the  method  is  to  obtain 
solutions  of  the  equations  of  motion  as  algebraic  expressions  with  unknown  constant  coefficients,  for  both  the 
slipping  and  the  stuck  states  of  the  joint.  Then,  to  calculate  those  coefficients,  two  types  of  conditions  are 
imposed:  a)  continuity  conditions  of  displacements  and  velocities  when  the  joint  passes  from  a  sliding 
situation  to  a  stuck  one,  and  b)  periodicity  conditions  relating  the  values  of  the  response  at  the  beginning  and 
at  the  end  of  a  semiperiod.  By  taking  advantage  of  the  form  of  the  resulting  equations,  the  mathematical 
problem  is  eventually  reduced  to  the  solution  of  a  system  of  algebraic  linear  equations. 

3.1  One  Sllp-to-Stuck  Change  In  a  Semiperiod. 

Without  any  loss  of  generality  we  can  define  the  initial  time,  t  =  0,  as  exactly  the  instant  at  which  the  joint 
passes  from  a  stuck  to  a  sliding  condition,  i.  e.,  the  instant  at  which  f  reaches  the  value  f0.  This  definition 
implies  the  introduction  of  an  unknown  phase  angle,  0,  in  the  terms  describing  the  excitation  forces  in  the 
equations  of  motion,  to  take  into  account  that  such  forces  are  not  known  at  this  initial  time.  Hence,  in 

equations  (5),  (6)  and  (8),  cos(cot)  and  cosc(cot  +  P)  must  be  replaced  by  cosjrot  +  <>)  and  cos(wt  +  p  +  0), 

respectively;  only  in  some  very  particular  cases  will  <j»  be  zero. 

The  tangential  force  f  while  the  joint  is  sliding  is  equal  to  the  constant  value  f0.  Under  this  condition 
equations  (6)  and  (7)  are  linear,  and  their  solutions  are 

x i (1)  =  Aifj(t)  +  Agfjft)  +  g,(t,0)  -  fg/(i-m)  (9) 

x2(t)  =  A3f3(t)  +  A4f4(t)  +  g2(t,<)>)  +  ff/m  (1 0) 

where  A;  (j=1 ,2,3,4)  are  unknown  constants,  fj(t)  (j=1 ,2,3,4)  are  known  functions  of  time  depending  on  the 
system  parameters  m,  k,  C,  and  C2,  and  gj(t,<>)  (j=1 ,2)  are  known  harmonic  functions  of  time  depending  on  the 
mentioned  parameters  and  those  defining  the  excitation  (a,  p,  and  a).  Functions  fj(t)  represent  the 
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complementary  solutions  of  the  equations,  and  functions  are  the  particular  solutions  due  the  harmonic 
excitation  All  of  them  are  given  in  Appendix  A. 

Differentiation  of  (9)  and  (10)  with  respect  to  time  yields  the  following  expressions  for  the  velocities: 
x,(t)  =  A,f,(t)  +  A2f2(t)  +  g,(t,<j>)  (11) 

*2(t)  =  Agf'ad)  +  A4f4(t)  +  g2(t,4>)  (12) 

The  first  condition  to  be  satisfied  by  this  solution  is  that  at  t  =  0  the  relative  velocity  between  the  two  bodies 
x^O)  -  x2(0)  is  zero,  in  agreement  with  the  definition  that  t  =  0  is  the  time  at  which  the  joint  is  passing  from  a 
stuck  to  a  sliding  condition.  Hence, 

A^fO)  +  A2f2(0)  -  A3f3(0)  -  A4f4(0)  =  g,(0,W  -  g2(0,4>)  (13) 

Simple  algebraic  operations  yield 

b, ,  A,  +  b12A2  +  bl3A3  +  b14A4  =  p1sin<>  +  q1cos<>  (14) 

The  coefficients  bjj,  p,  and  q,  are  defined  in  Appendix  A.  Note  that  these  are  time  independent  quantities. 

We  will  now  assume  that  the  joint  remains  sliding  during  some  time  interval  tt  >  0  at  the  end  of  which  the 
relative  velocity  becomes  zero  again  and  the  joint  returns  to  a  stuck  state.  This  is  expressed  by 

A,it(t,)+  A2f2(t,)  -  A3f3(t,)  -  A4f4(tf)  =  g,(t(,0)  -  g2(t(,4>)  (15) 

that  can  be  also  written  in  the  form: 

b2j  A^  +  b22A2  +  b23A3  +  b24A4  —  p2sintj)  +  q2cos0  (15) 

In  addition,  we  assume  that  between  t,  and  the  end  of  the  semiperiod  t2  (t2  *  rt/m)  the  joint  remains  stuck. 
During  this  second  time  interval  the  motion  of  the  system  is  described  by  the  linear  ordinary  differential 
equation  (8),  whose  solution  is 

x,(t)  =  E^sd)  +  E2f6(t)  +  g3(t,<t>)  +  kd  (17) 

As  before,  f5(t)  and  f6(t)  are  known  complementary  solutions  of  the  differential  equation  which  depend  on  the 
system  parameters,  and  g3(t,o)  is  the  particular  solution  due  to  the  harmonic  external  forces  which  is 
determined  by  both  the  system  and  the  excitation  parameters.  These  functions  are  also  described  in 
Appendix  A. 

The  expression  for  the  common  velocity  of  the  two  bodies  is  now 

x1(t)  =  E1f5(t)  +  E2f6(t)+g3(tD)  (18) 

The  constants  d,  E,  and  F2  can  be  found  from  the  continuity  conditions  of  the  displacements  and  velocities  at 
t  =  t, ,  which  require  these  quantities  to  have  respectively  the  same  values  just  before  and  just  after  t,  In  the 
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case  of  x,(t),  replacing  I  by  t,  ;nlo  (9)  and  (17)  and  equating  the  results  we  obtain 
E,f5(t,)  +  Ejfg(t^)  +  g3(t,.<W  +  kd  =  A1f  1  (t-| )  +  AgfgOi )  +  Q 1 0 1  >*!*)  ‘  ty(1-m) 
Similarly,  for  x2(t),  from  (7),  (12)  and  (17)  we  have 

Ei^t,)  +  Ejfg^,)  +  gg(t^  ,<}>)  +  kd-d  =  Agfg(t^)  +  A4f4(t1 )  +  g2(t-|,<}>)  + 
and  for  x^t), 

^1^5(1 1 )  +  +  93(^1  *4*)  =  A2f  1  (t  1 )  +  Ag^tj)  +  3i  (*1  -4*) 


The  continuity  condition  for  x2(t)  is  automatically  satisfied  uy  imposing  it  to  x^t),  since  t,  is,  by  definition,  a 
time  at  which  both  velocities  have  the  same  value. 

The  last  three  equations  constitute  a  linear  system  in  d,  E1  and  E2,  whose  solution  (Appendix  B)  is 
d  =  A^f ^ (t ^ )  +  A2f2(t1  +  9](t^,<t))  -  -m)  -  Agfg(t^)  -  A4f4(t1)  -  g2(t.|  >4)  '  (22) 

E,  =  A1f7(t1)  +  A2fg(t1)  +  Agfg(t,)  +  A4f10(t()  +  g4(t  !•<())  +  C^q  (23) 

E2  =  A1f11(t1)  +  A2f12(t1)  +  Agf^gjt^  +  A4fj4(tj)  +  95(^1 .4*)  +  C2(^i^o  (24) 


As  the  excitation  is  harmonic,  the  steady-state  response  at  a  time  t  has  to  have  the  same  magnitude,  but 
opposite  sign,  than  at  t  plus  or  minus  a  semiperiod,  t2  .  This  requirement  is  referred  as  the  periodicity 
condition  and  can  be  expressed  as 

x1(t2)  =  -x1(0)  (25) 

x2(t2)  =  -  x2(0)  (26) 

x1(t2)=-x1(0)  (27) 

f(t2)  -  - 1(0)  (28) 

The  last  condition  is  feasible  because  at  t  =  0,  and  as  a  consequence  at  t  =  t2,  the  force  in  the  joint  is  a 
.  ntinuous  function  of  time.  Later  on,  we  will  see  that  this  is  not  the  case  at  t  =  t,.  The  periodicity  condition 
for  x2(t)  coincides  with  the  one  for  x,(t)  because  both  velocities  are  identical  at  t=0  and  also  at  t  - 12. 

Equation  (25)  can  be  expanded  as 

A-t f  1  (0)  +  A2f?(0)  +  g^(0,(t>)  -  fg/(1-m)  +  E^fg(t2)  +  E2fg(t2)  +  93(^2’*^)  +  kd  =  0 
Using  the  expressions  previously  obtained  for  d,  E,,  and  E2  we  arrive  at 

bg^A^  +  bggAg  4-  bggAg  +  bg4A4  +  bggfg  =  PgSih({)  +  qgCOS^I  (29) 

Similarly,  (26)  and  (27)  can  be  written  as 
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bjiA,  +  b42A2  +  b43A3  +  b44A4  +  b43fQ  =  p4sin$  +  q4cos<t> 
b5iA,  +  b52A2  +  b53A3  +  b54A4  +  b5510  =  pssirvt>  +  q5cos$ 


(30) 

(31) 


The  last  peridiocity  condition,  equation  (28),  says  that  at  the  end  of  the  semiperiod  the  tangential  force  in 
the  joint  is  equal  to  -f0  .  In  fact,  at  exactly  this  time  another  change  from  the  stuck  to  the  sliding  condition 
occurs,  as  it  happens  at  t  =  0.  The  mentioned  force  is  computed  by  replacing  xt(t)  from  (17)  into  the  proper 
terms  of  the  differential  equations  (5)  or  (6),  with  t  *  t2.  The  result  can  be  written  as: 

b6 1 A  i  +  bgjAj  +  b63A3  +  p6sind  +  qgCOSit.  =  0  (32) 

In  summary,  we  have  formulated  the  conditions  to  find  an  analytical  solution  for  the  steady-state  vibrations  of 
the  system  under  study  when  only  one  slip-to-stuck  change  occurs  in  the  friction  joint  in  a  semiperiod  of 
vibration.  The  result  is  a  set  of  six  coupled  algebraic  equations  ((14),  (16)  and  (29)  to  (32))  with  six  unknowns 
(A,  to  A4,  t,  and  4>).  Although  it  is  mathematically  possible  to  solve  directly  this  system  of  equations,  this  is  a 
complicated  task  because  most  of  the  coefficients  bjj  are  strongly  nonlinear  functions  of  the  unknown  t,  and 
some  iterative  numerical  procedure  would  have  to  be  used. 

An  indirect,  simpler  approach  is  used  here  taking  into  account  that  the  first  five  equations  of  the  system  are 
linear  in  f0.  We  consider  f0  as  one  of  the  unknowns  and,  in  exchange,  tt  as  an  independent  variable.  In  this 
way,  the  coefficients  b,(  become  constants  and  the  first  five  equations  can  be  written  as 

B  A  =  Psind  +  Qcosd  (33) 

where 

A  —  A2  A3  A4  fg  > 

b,2  ...  b15 

l 

b22  ...  b25 

i 

i 

b52  -  b55 

PT  =  <p,  p2  p3  p4  p5  > 

Qt  =  <q,  q2  q3  q4  q5  > 

Solving  this  linear  system  of  equations  we  obtain 

A  *  B_1Psind  +  B^Qcosd  (34) 
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in  particular, 

A,  =  p7sir»4>  +  q^sa 
a2  =  p8sin(j)  +  qgcoaj) 

A3  =  pgsin<j)  +  qgCOS<!> 

Replacing  these  values  into  (32),  we  arrive  at 

(k>6i  P7  +  b62P8  +  b63P9  +  Ps^®*0^  +  (^61^7  +  ^62%  +  b63P9  +  P^COSij)  =  ® 
or 

,  ,  b6lP7  +  b62Ps  +  b63P9  +  P6 

-tan<j>  =  - - - - - - 

b6lP7  +  b62P8  +  b63*^9  +  ^6 


(35) 


This  equation  yields  two  values  of  4>.  For  each  of  them  a  set  of  values  for  the  elements  of  vector  A  is 
obtained  with  (33).  Each  of  these  sets  is  the  negative  of  the  other,  reflecting  the  periodicity  condition.  We  are 
interested  in  the  solution  with  a  positive  f0,  because  we  have  assumed  that  this  is  the  case  at  t=0.  The  other 
solution  is  equally  valid  and  corresponds  to  the  other  semiperiod  of  vibration. 

It  should  be  pointed  out  that  the  above  solution  is  valid  only  if  the  relative  velocity  x^t)  -  x2(t)  does  not  have 
any  change  of  sign  between  0  and  f,;  otherwise,  there  will  be  several  changes  from  the  sliding  to  the  stuck 
condition  and  viceversa,  contradicting  the  assumption  of  only  one  of  such  changes  in  a  semiperiod.  This 
requirement  was  not  imposed  when  formulating  the  equations;  therefore,  it  has  to  be  verified  separately.  For 
the  same  reasons,  it  also  has  to  be  verified  that  the  magnitude  of  the  tangential  force  in  the  friction  joint 
between  tt  and  t2  is  smaller  than  f0. 

3.2  Friction  Joint  Always  Sliding 

When  the  normal  load  is  smaller  than  some  limit  value  the  friction  joint  is  always  sliding.  In  this  case  t,  =  t2 
and  the  differential  equations  (5)  and  (6),  as  well  as  their  solutions  (9)  and  (10),  are  valid  all  the  time. 
Assuming  that  a  change  of  sign  in  the  relative  velocity  and,  as  a  consequence,  in  the  friction  force,  occurs  at 
t=0,  which  requires  the  introduction  of  the  unknown  phase  angle  <t>  in  the  excitation  force,  we  have  now  five 
unknowns,  namely  A,,  A2,  Ag,  A4  and  $.  This  assumption  leads  again  to  equation  (14). 

The  periodicity  conditions  are  also  expressed  by  equations  (25)  to  (27)  and  by 
x2(t2)  =  -  x2(0)  (36) 

which  this  time  is  not  automatically  satisfied  because  we  have  not  enforced  the  condition  that  the  velocities 
are  the  same  at  t  =  t2  (this  is  an  equivalent  requirement). 
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Equations  (25)  and  (27)  can  be  written  again  in  the  form  of  (29)  and  (31).  However,  an  additional 
simplification  can  be  done  considering  that  gj(t,$)  are  harmonic  functions  with  period  2t2  and,  as  a 
consequence,  gj(t2,4>)  =  -gj(0,<(>).  The  same  is  true  for  any  time  derivatives  of  gj(t,<M-  Therefore,  instead  of  (29) 
«md  (31),  we  now  have 

^,(0)  +  AgfjtO)  -  2f(j/(1-m)  +  A,f ,(t2)  +  A2f2(t2)  -  0  (37) 

A^tO)  +  A2f2(0)  +  A1f1(t2)  +  A2f2(t2)  =  0  (38) 

This  constitutes  a  system  of  two  linear  equations  with  two  unknowns,  A,  and  A2,  that  can  be  solved 
explicitely. 

Following  the  same  procedure  with  (26)  and  (36),  a  similar  system  of  equations  for  A3  and  A4  is  obtained. 
After  all  the  Aj  are  known  <j>  can  be  calculated  with  (14).  As  in  the  previous  case,  the  condition  that  no  change 
in  the  relative  velocity  occurs  between  t=0  and  t=*t2  has  to  be  verified  separately. 

3.3  Some  Limiting  Cases 

The  analytical  solutions  presented  in  this  section  provide  the  tools  for  computing  the  values  of  the 
parameter  f0  corresponding  to  three  limiting  situations.  First,  recalling  that  t,  was  defined  as  the  time  at  which 
slip  occurs  in  a  semiperiod  of  vibration,  when  t,«0.  there  is  no  slip  at  all,  and  equation  (8)  with  d  =  0  always 
governs  the  motion  of  the  system.  Therefore,  the  solution  for  this  case,  gives  for  pN/Q,  a  limit  r,  such  that 
when  uN/Q,  >  r,  the  joint  remains  always  stuck. 

The  second  limiting  case  is  that  of  no  normal  force  in  the  joint,  i.e.,  f0  «  0  .  The  joint  is  always  sliding  and 
the  two  masses  vibrate  separately.  The  solution  can  be  obtained  introducing  f0  =  0  in  equations  (9)  and  (10). 

Finally,  we  can  obtain  the  maximun  value,  r2,  of  pN/Q1  for  which  the  joint  will  be  always  sliding.  This  limit 
can  be  calculated  by  considering  t1  equal  to  the  semiperiod  of  vibration  (t1  »  n/ca)  in  the  equations  of  section 
3.1.  The  same  result  is  obtained  solving  the  equations  for  the  case  of  pure  slip  (section  3.2)  with  the 
additional  restriction  that  the  initial  value  of  the  relative  acceleration  is  zero.  The  reason  is  that  our  solutions 
are  valid  only  if  the  relative  velocity  is  positive  during  the  first  semiperiod  of  vibration,  and  since  the  initial 
velocity  is  nil,  zero  is  the  lower  value  of  the  relative  acceleration  at  t  =  0  that  can  satisfy  this  condition. 
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3.4  Numerical  Integration 

Theoretically,  the  analytical  method  of  solution  can  be  easily  extended  to  cases  with  more  slip-to-stuck  and 
stuck-to-slip  changes  in  the  friction  interface,  by  using  in  each  situation  the  analytical  solution  of  the 
corresponding  differential  equations,  and  by  imposing  similar  continuity  and  periodicity  conditions.  However, 
the  size  and  complexity  of  the  resulting  algebraic  problem  increase  very  rapidly,  and  numerical  step-by-step 
integration  of  the  differential  equations  becomes  a  more  practical  method  to  obtain  an  exact  solution. 

Several  well  known  time  step-by-step  schemes  can  be  used  to  integrate  the  equations  of  motion,  for 
instance,  the  Runge-Kuta  or  the  Newmark  families  of  methods.  Efficient  algorithms  to  apply  these  methods, 
as  well  as  discussions  of  their  mathematical  characteristics  (stability,  convergency  and  accuracy),  are 
presented  in  a  number  of  books  on  numerical  methods  and  on  structural  dynamics  [11,12],  More  recently, 
Dickens  and  Wilson  [13],  Sinha  and  Griffin  [4],  Hilber  et  al  [14],  Bazzi  and  Anderheggen  [15],  and  Zienkiewicz 
et  al  [16]  have  presented  improved  versions  of  these  type  of  methods. 

For  the  problem  considered  in  this  paper  almost  all  the  available  methods  could  be  used  with  more  or  less 
the  same  levels  of  efficiency  and  accuracy,  because  we  are  dealing  with  only  two  ordinary  differential 
equations.  Obviously,  this  remark  is  not  true  tor  problems  involving  more  degrees  of  freedom  and  more 
friction  joints. 

The  main  advantage  of  using  step-by-step  procedures  is  that  any  number  of  slip-to-stick  changes  and  vice 
versa  can  be  considered  without  difficulty.  In  fact  that  number  is  a  result  of  the  solution  procedure.  On  the 
other  hand,  the  main  disadvantage  arises  from  the  fact  that  for  nonlinear  steady-state  problems  the 
integration  has  to  be  carried  out  for  a  sufficiently  long  time  to  reach  the  periodic  state,  i.e.,  until  the  effect  of 
the  initial  disturbances  becomes  negligible.  For  this  reason  the  numerical  integration  can  be  impractical  even 
for  systems  with  a  moderate  number  of  degrees  of  freedom.  In  this  work  this  approach  was  used  mainly  to 
verify  the  analytical  exact  method  and  to  assess  the  accuracy  of  approximate  methods  in  cases  with  several 
slip-stick  changes  in  the  friction  interface  (Fig.2). 

3.5  Numerical  Results 

In  order  to  understand  the  role  of  the  different  parameters  defining  the  system  depicted  in  Fig.  i  in  its 
dynamic  behavior,  we  have  computed  the  response  for  particular  values  of  those  parameters.  Of  special 
interest  is  the  influence  of  the  friction  joint  in  the  response.  We  also  want  to  gain  insight  about  the  usefulness 
of  analytical  results  in  designing  structures  with  friction  interfaces. 

The  selected  values  of  the  parameters  are  presented  in  Table  1 .  Detailed  exact  results  for  the  first  case  of 
this  table  are  presented  in  this  section.  We  have  verified  that  all  the  relevant  characteristics  of  the  steady- 
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state  response  of  th8  system  under  study  can  be  appreciated  in  the  results  corresponding  to  this  case. 
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me  vaiues  of  the  ratio,  rv  of  the  maximum  frictional  force,  pN,  to  the  amplitude,  Qv  of  the  the  load  applied 
to  M,  that  prevent  slipping  are  presented  in  Fig.  3  as  a  function  of  the  excitation  frequency.  For  a  prescribed 
frequency,  any  pN/Q,  ratio  equal  or  greater  than  the  value  in  the  curve  will  not  allow  any  slipping  in  the  joint. 
On  the  other  hand,  a  prescribed  value  of  pN/Qv  can  be  represented  by  a  horizontal  line  whose  intersections 
with  the  curve  of  limiting  values  provide  the  range  of  frequencies  for  which  there  will  be  some  slipping  in  the 
joint  The  peak  value  (12.96)  of  this  curve  occurs  for  the  resonant  frequency  of  the  system  when  the  joint 
remains  locked  (o=1).  The  maximum  value,  r2,  of  pN/Q,  that  allows  the  joint  to  be  sliding  all  the  time  is  also 
presented  in  Fig.  3.  For  a  given  frequency,  the  joint  will  be  always  sliding  if  pN/Q,  is  equal  or  lower  than  the 
limiting  value  in  the  curve;  and  for  a  prescribed  pN/Q,  we  can  determine  the  range  of  frequencies  for  which 
the  joint  is  always  slipping. 

The  effects  of  pN/Q,  on  the  amplitude  of  the  system  response  are  illustrated  in  Figs.  4,  5  and  6,  in  which 
jvquency  response  curves  have  been  plotted  for  several  values  of  pN/Q,  .  Fig.  4  shows  that  the  amplitude, 
_.  ,x,  of  the  normalized  displacement  of  M<,  has  two  sets  of  peak  values.  The  first  one  corresponds  to  the 
normalized  natural  frequency  of  the  Mj-K,  system,  on  which  the  external  load  is  applied.  The  maximum, 
resonant,  response  occurs  when  there  is  no  friction  force  (pN/Q^O)  and  M1-K1  is  vibrating  without  any 
interaction  with  the  other  mass  and  spring.  For  small  values  of  pN/Q,  the  friction  joint  is  always  slipping  and 
ihe  peak  values  occur  for  the  same  resonant  frequency.  However,  the  presence  of  these  small  friction  forces 
i ..  rduces  appreciable  reductions  of  the  maximum  response. 

As  pN/Q,  increases,  the  peak  values  move  to  an  excitation  frequency  close  to  unity,  the  dimensionless 
resonance  frequency  of  the  system  when  the  joint  remains  locked.  In  such  cases,  slipping  occurs  during  one 
part  of  a  cycle  only.  Eventually,  if  pN/Q,  equals  or  exceeds  the  limiting  value  presented  in  Fig.  3  (12.96  for 
m=l),  the  two  masses  will  vibrate  together  without  any  relative  slipping  (compare  Figs.  5  and  6  for 
pN/Q, =12.96).  This  process  can  be  appreciated  by  examining  the  detailed  time  histories  depicted  in  Figs.  7 
and  8.  In  Fig.  7,  pN/Q,  is  large  enough  to  prevent  slipping.  Thus  x,  and  x2  coincide,  and  both  have  the 
maximum  value  51 .7.  This  is  an  elastic  case  and  all  the  steady-state  responses  are  harmonic.  In  Fig.  8,  with 
a  lower  value  of  pN/Q,  ,  there  is  a  small  difference  between  x1  and  x2,  resulting  from  a  small  slip  between  the 
two  masses  (shown  by  a  dashed  line).  This  slip,  however,  is  sufficient  to  produce  a  significant  reduction  in  the 
maximum  values,  which  are  now  30.5  and  31 .3,  respectively.  Even  though  in  this  case  the  displacements  are 
not  strictly  harmonic,  it  is  obvious  that  they  can  be  accurately  approximated  by  harmonic  functions  with  the 
same  frequency  as  the  excitation. 

In  the  case  of  Fig.  9  the  value  of  pN/Q,  is  very  small  and  consequently  the  joint  is  always  sliding.  The 
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friction  force  varies  as  an  alternating  square  wave.  As  in  the  previous  case,  the  variations  of  the 
displacements  are  esentially,  although  not  exactly,  harmonic  waves. 

The  main  purpose  of  including  a  friction  element  in  a  vibrating  system  is  to  dissipate  some  of  the  energy 
imparted  to  the  system.  Since  energy  dissipation  is  a  function  of  the  slip  load  pN  and  the  fraction  of  the 
period  of  oscillation  during  which  slip  occurs  it  is  of  interest  to  examine  how  this  fraction,  denoted  as  x  in  Fig 
1 0.  varies  with  the  frequency  of  excitation  co  and  the  normalized  slip  load  pN/Q,.  x  is  shown  in  Fig.  1 1 .  Note 
that  for  most  values  of  pN/Q,  and  a  the  friction  element  either  slips  during  the  entire  period  of  oscillation  (x=l) 
or  remains  stuck  (x=0).  The  transition  generally  occurs  over  a  narrow  frequency  band.  In  particular,  for  small 
values  of  pN/Q1  and  frequencies  greater  than  unity  a  sudden  change  from  a  totally  sliding  condition  to  one  of 
no  sliD  can  be  appreciated.  This  effect  may  be  of  significance  in  design. 

For  cases  of  pure  slip,  as  that  of  Fig.  9,  the  tangential  force  in  the  friction  element  changes  abruptly  from 
+pN  to  -pN  as  the  relative  velocity  between  the  end  points  of  the  element  changes  sign.  Clearly,  no  jump 
occurs  if  the  element  is  stuck  (Fig.  7).  For  partial  slip  the  jump  takes  a  value  that  lies  between  those 
corresponding  to  the  extreme  cases  (Fig.  8).  The  values,  q,  of  this  jump  in  the  slip  force,  normalized  with 
respect  to  pN  (see  Fig.  1 1)  are  presented  in  Fig.  12.  Note  the  similarity  between  the  x  and  q  curves. 

It  can  be  shown  that  the  energy  dissipated  in  each  cycle  of  vibration  is  equal  to  4  times  the  product  of 
uN/Q,  and  the  initial  relative  displacement  x,(0)-x2(0).  Obviously,  the  energy  is  zero  when  either  the  relative 
displacement  or  the  normal  force  in  the  joint  are  nil,  and  it  will  have  a  maximum  value  for  some  intermediate 
value  of  pN/Q,  .  This  can  be  appreciated  in  Fig.  13  for  (o=0.7071  and  in  Fig.  14  for  to=1.  The  variation  of 
energy  with  the  excitation  frequency,  for  prescribed  values  of  pN/Q,  ,  is  shown  in  Fig.  15. 

The  phase  angles  y,  and  y2,  of  the  displacements  x1  and  x2  with  respect  to  the  exciting  force  are  shown  in 
figures  16  and  17  as  a  fraction  of  n.  Two  elastic  limiting  cases  are  included,  namely  the  case  of  no  friction 
(pN/Q,=0)  and  that  of  no  slip  (pN/Q,  =12.96).  Fig.  16  shows  that  y,  is  always  positive,  indicating  that  the 
response  always  lags  behind  the  excitation,  as  in  the  elastic  case.  For  frequencies  between  the  resonant 
ones  (i.e.,  between  the  almost  vertical  slopes  of  the  limiting  cases)  y,  lies  between  the  two  elastic  values 
corresponding  to  that  frequency.  For  y2  we  have  only  one  limiting  elastic  case  because  there  is  no  force 
acting  on  mass  M2. 

Fig.  18  shows  the  maximum  values,  x1max  ,  of  the  normalized  displacement  of  mass  M,,  as  a  function  of 
pN/Q,.  All  the  possible  values  of  the  excitation  frequency  were  considered  in  obtaining  each  point  of  this 
curve.  For  low  values  of  pN/Q,  ,  x1max  decreases  very  rapidly,  until  a  minimum  value  is  reached  (x1max  = 
6.36)  Then,  for  further  increments  of  pN/Q,,  x1fnax  also  increases,  until  a  constant  value  is  obtained  when 


12 


>r,Vw\ 


pN/Q,  is  large  enough  to  mantain  the  joint  always  locked.  The  value  of  pN/Q,  corresponding  to  the  minimum 
response  is  defined  as  the  optimum  value  of  such  a  ratio  (in  this  case,  pN/Q,  optimum  =  1.0).  Similar 
optimization  curves  can  be  constructed  for  any  other  response  of  interest.  Naturally,  it  is  highly  desirable  to 
design  the  friction  joint  to  have  parameters  close  to  the  optimum  ones. 

In  order  to  evaluate  the  accuracy  that  can  be  expected  from  approximate  methods,  we  have  computed  the 
coefficient  of  the  first  Fourier  component  of  the  tangential  force  in  the  joint.  Results  for  co=0.707l  and  m=l  are 
presented  in  Figs.  19  and  20.  As  expected,  these  coefficients  have  a  maximum  value  of  4 In  corresponding  to 
cases  in  which  the  joint  is  always  sliding.  The  minimum  value  is  1 .0  and  occurs  when  the  joint  is  locked. 

Exact  results  have  been  also  computed  for  other  cases  included  in  Table  1.  As  pointed  out  earlier,  they 
present  characteristics  similar  to  those  discussed  in  the  case  covered  in  this  section.  These  exact  results  are 
later  used  in  comparisons  with  results  of  approximate  methods. 

4  APPROXIMATE  METHODS  OF  SOLUTION 

In  view  ot  the  difficulties  of  obtaining  analytical  or  numerical  exact  solutions,  approximate  methods  have 
been  used  by  several  researchers  [4,  6,  7],  The  most  widely  used  method  is  the  so-called  harmonic  balance 
method  which  finds  the  best  approximation  of  the  response  in  the  space  of  the  harmonic  functions  with  the 
same  frequency  as  the  excitation.  This  method  is  completely  equivalent  to  the  linerization  technique 
proposed  by  Caughey  [17],  Its  theoretical  foundations  have  been  thoroughly  studied  by  Spanos  and  Iwan 
[18]. 

in  this  section  a  complex  variable  version  of  the  method  is  presented  for  the  system  under  study;  then, 
some  modifications  are  proposed  to  improve  its  accuracy  when  the  joint  is  not  always  sliding.  A  comparison 
of  exact  and  approximate  results  is  also  presented  to  ilustrate  the  accuracy  and  applicability  of  the  various 
methods. 

4.1  Harmonic  Balance  Method 

The  following  assumptions  are  made  in  this  method: 

a)  The  response  is  harmonic  with  the  same  frequency  as  the  excitation  forces; 

b)  The  friction  force  has  the  same  direction  as  the  relative  velocity  in  the  joint,  but  opposite  sign; 

c)  The  friction  joint  is  always  sliding,  i.e.,  no  cuck  condition  is  considered; 

d)  Only  the  first  Fourier  component  of  the  friction  force  has  a  significant  participation  in  the 
response  of  the  system,  i.e.,  the  effect  of  higher  components  is  negligible. 

In  order  to  present  a  complex  variable  app'oach  of  the  method  the  excitation  forces  are  defined  as 
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Q/t)  =  qexp(iWT) ,  j=1.2 

where  an  overbar  denotes  a  complex  quantity,  and  i=V=T 
According  to  assumptions  a)  and  d)  we  can  also  write 


Xj(T)  =  Xjexp(iWT)  (40) 

F(T)  =  F  exp(iWT)  (41) 

thus, 

Xj(T)  =  iW  Xj  exp(iWT)  (42) 

Xj(T)  =  -  W2Xj  exp(iWT)  (43) 


By  substituting  these  expressions  into  (1)  and  (2)  and  dropping  the  factoi  exp(iWT)  the  following  complex 


algebraic  equations  are  obtained: 

(Kr  W2M,  +  iWC,)  X,  =  Q,  -  F  (44) 

(K2-  W2M2  +  iWC2)  X2  .  Q2  +  F  (45) 

solving  for  Xt  and  X2 

X,  =  ^(WJQ,  -  h,(W)F  (46) 

x2  =  h>m  Q2  +  tyW)  F  (47) 

where 

f)(W)  =  (Kj  -  W2Mj  +  iWCj)'1  (48) 

The  relative  displacement  can  be  now  written  as 

X,  -  X2  =  (h,(W)  Qr  h2( W)  Q2)  -  (h,(W)  +  h,( W)JF  (49) 

define 

D  =  D  exp  (i0)  =  X,-  X2  (50) 

P=  ft(W)Q,  -  h,(W)  Q2j  (51) 

R  =  (h,(W)  +  fk(W)]  (52) 


where  D  is  the  amplitude  and  0  is  the  phase  angle  of  the  relative  displacement  in  the  friction  joint.  This  allows 
us  to  write 

D  +  RF  =  P  (53) 

which  is  an  equation  relating  the  relative  displacement  with  the  tangential  force  in  the  friction  joint. 

From  (42),  the  relative  velocity  is 

X,(T)  -  X2(T)  =  iW(Xr  X2)exp(iWT)  =  iWD  exp(iWT) 

Taking  into  account  that  W  is  a  positive  real  constant,  the  definition  of  0  in  (50),  and  that  i  =  exp(jt/2),  we 
conclude  that  the  phase  angle  of  the  relative  velocity  in  the  joint  is  (0  +  «/2).  According  to  assumption  b)  this 
is  also  the  phase  angle  of  the  friction  force  (the  proper  sign  was  already  considered  when  formulating  the 
equations  of  motion).  Therefore, 

F  =  iF  exp(i0)  (54) 


Besides,  according  to  assumptions  c)  and  d)  this  force  has  the  shape  of  an  alternating  square  wave  in  time 

(see  Fig.  4),  and  only  the  first  Fourier  coefficient  of  such  a  wave  has  to  be  considered.  Hence, 

F  =  (4/jt)pN 

Then,  for  a  prescribed  value  of  pN,  the  only  unknows  in  (53)  are  D  and  0. 

The  traditional  application  of  the  harmonic  balance  method  consits  in  equating  the  real  and  the  imaginary 
parts  of  both  members  of  (53)  obtaining  two  real  nonlinear  equations  in  D  and  0.  The  equation  resulting  from 
the  real  parts  corresponds  to  a  cosine  excitation  force,  and  the  one  resulting  from  the  imaginary  parts 
corresponds  to  a  sine  excitation.  Any  of  them  can  be  expressed  as 
c(W,D,F,0)  cos(WT)  +  s(W,D,F,0)  sinfWT)  =  0. 

Since  this  equation  has  to  be  satisfied  for  any  value  of  T,  both  c(W,D,F,0)  and  s(W,D,F,0)  have  to  be 
separately  equal  to  zero,  providing  two  equations  for  D  and  0.  Elimination  of  0  leads  to  the  following  quadratic 
equation  in  D: 

a,  (W)D  +  a2(W)DpN  +  a3(W)  (pN)2  =  a4(W)  (55) 


where  aj(W)  (j=1,4)  are  frequency  dependent  coefficients. 

After  solving  for  D,  0  is  calculated  with  one  of  the  two  original  equations.  Then,  (54),  (46)  and  (47)  are  used 


to  compute  F,  X,  and  X2. 

An  easier  way  to  obtain  (55)  is  to  multiply  (53)  by  it  complex  conjugate;  the  result  is 

DD*  +  DRT  +  D*RF  +  RRTF*  =  PP*  (56) 

where  the  star  denotes  the  complex  conjugate. 

From  (50)  and  (54)  we  have: 

D*  =  D  exp(-iG) 

T  =  -iF  exp(-iG) 
and  consequently 

DF'  =  -iDF  exp(i0)  exp(-iG)  =  -  iDF 

dT  =  iDF  exp(iG)  exp(-iG)  =  iDF 

DR*F*  +  o¥  =  -iDFR*  +  iDFR  =  iDF[R-  R*j  =  iDF[2i][Imag(R)j=  -2DF[lmag(R)] 

(57) 

where  Imag(R)  denotes  the  imaginary  part  of  R.  Recalling  that  the  product  of  a  complex  quantity  by  its 
conjugate  yields  the  square  of  its  modulus  we  also  have  DD*  =  D2,  RR*  =  R2,  FF*  =  F2,  and  PF5  =  P2. 
Therefore,  (56)  becomes 

D2  -  2  Imag(R)  FD  +  R2F2  =  P2  (58) 

which,  as  (55),  is  a  quadratic  equation  for  D. 

4.2  Modified  Harmonic  Balance  Methods 

As  it  could  be  expected,  the  harmonic  balance  method  yields  very  accurate  results  when  the  variation  of 
the  friction  force  with  time  is  similar  to  the  alternating  square  wave  shown  in  Fig.  9.  Larger  errors  occur  when 
the  joint  is  always  locked  (Fig.  7),  because  the  force  in  the  joint  is  actually  harmonic  and  the  value  of  its  first 
and  only  nonzero  Fourier  coefficient  is  1.0  and  not  4/ji,  or  when  slipping  is  prevented  most  of  the  time  (Fig.  8), 
because  such  a  coefficient  will  be  closer  to  1  than  to  4 In.  Appreciable  errors  also  occur  when  one  or  both  of 
the  displacements  have  significant  components  of  higher  frequencies,  as  in  the  example  shown  in  Fig.  23. 

Aiming  to  reduce  the  first  of  the  above  sources  of  error,  some  modifications  to  the  method  are  proposed 
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here.  They  consist  in  estimating  the  first  Fourier  coefficient  of  the  friction  force  with  heuristic  rules,  instead  of 
considering  it  always  equal  to  4 In.  Thus  the  amplitude  of  F,  is  now  given  by 

F  =  C,pN  (59) 

where  c,  is  estimated  with  some  rule  satisfying  the  condition 

1 .0  <  C,  <  4/rt  (60) 

The  rules  proposed  here  have  the  form 

Cf=  C,(D)  (61) 

which  means  that  the  coefficient  is  considered  to  be  a  function  of  the  amplitude,  D,  of  the  relative 
displacement.  Obviously,  Cf(D)  is  1 .0  when  D  =  0,  and  its  maximum  value  is  4/?c . 

The  first  rule  is  expressed  in  terms  of  a  dimensionless  parameter  0,  and  is  called  parametric.  C,  is  given  by 
C.(0)  =  £2!  +  2sin0]  ht  (62) 

1  COSv 

o  is  related  to  D  by 

D  =  D1(29  +  sin  20)/jc  (63) 

where  D1  is  obtained  by  substituting  F  =  pN  into  (58).  The  variation  of  C,(0)  with  D,  according  to  this  rule  is 
shown  in  Fig  21. 

The  next  step  is  to  replace  (59),  (62)  and  (63)  into  (58)  to  obtain  a  nonlinear  equation  for  0.  After  solving 
numerically  this  equation,  (59)  and  (63)  are  used  to  calculate  D  and  F. 

A  second,  simpler  rule,  also  illustrated  in  Fig.  21,  consists  in  assuming  that  C,  is  a  linear  function  of  D,  with 
values  between  1.0  for  D  =  0  and  4/jt  =  1.2732  for  D  =  D1,  i.e.: 

Cf  =  1 .  +  0.2732  (D/D1)  (64) 


This  rule  has  the  advantage  that  substitution  of  F  by  C,pN  into  (58)  yields  another  quadratic  equation,  without 
adding  complexity  to  the  algebraic  problem.  This  approach  is  referred  as  the  modified  linear  method. 

A  third  rule,  similar  to  the  second  one,  consists  in  prescribing  a  parabolic  variation  for  Cf  within  the  interval 
0  and  D1,  as  it  is  also  shown  in  Fig.  21 .  In  this  case 

Cf  =  1. +0.2732  (D/D1)2  (65) 

Substitution  into  (58)  produces  this  time  a  cubic  real  equation  in  D.  This  option  will  be  called  modified 
quadratic  method. 
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A  disadvantage  of  the  preceeding  rules  is  that  equation  (58)  has  to  be  solved  twice  for  a  given  value  of  pN, 
first  with  C,  =  1 ,  to  compute  D1,  and  then,  for  the  corresponding  estimate  value  of  Cf.  In  the  present  study  the 
additional  numerical  work  is  not  excessive;  but  it  can  be  very  appreciable  in  problems  with  multiple  joints. 

Reductions  in  the  numerical  work  when  using  the  second  and  third  rules  can  be  obtained  if  D/D1  is 
computed  assuming  that  the  variation  of  D  with  F  is  linear,  as  shown  in  Fig.  22.  First,  the  value  of  pN 
corresponding  to  D  =  0  is  calculated  taking  into  account  that  in  this  case  the  Fourier  coefficient  of  the  force  in 
the  friction  joint  is  1 .  From  (58),  this  value  is: 

F  =  pN0  =  P/R  (66) 

For  any  value  of  the  amplitude  of  the  friction  force  equal  or  greater  than  pN0  the  joint  will  be  locked  and  the 
solution  will  remain  the  same.  This  provides  an  exact  criterion  to  define  the  slip-stick  limit  for  the  joint 
because  when  the  joint  is  locked  the  response  is  actually  harmonic.  It  is  important  to  take  into  account  that  it 
is  possible  to  obtain  physically  acceptable  solutions  only  for  values  of  pN  lower  than  pN0.  Now,  a  linear 
variation  of  D  with  pN  can  be  expressed  as 

D  =  D0(1-pN/pN0)  (67) 

.vnere  D0  is  the  amplitude  of  the  relative  displacement  corresponding  to  no  friction  force  in  the  joint. 

For  the  modified  linear  method,  combinig  (64)  and  (67)  we  obtain 

D/D1  =  1/(1  +  2732y)  (68) 

where  y  is  a  nondimensional  parameter  defined  as 

y  =  (pN/pN0)/[  1  -  pN'pNJ  (69) 

Therefore,  we  can  compute  C,  directly  and  then  solve  equation  (58).  Note  that  it  is  not  necessary  to  compute 
the  value  of  D0. 

For  the  quadratic  modification,  from  (65)  and  (67)  we  obtain 

D/D1  =  (1  +  1 ,093y  -  1  )/(.5465y)  (70) 

and  we  can  proceed  as  in  the  previous  case. 

Again,  the  proposed  modifications  are  intended  to  improve  the  perfomance  of  the  harmonic  balance 
method  when  higher  frequencies  do  not  play  a  significant  role  in  the  system  response.  If  this  is  not  the  case, 
appreciable  differences  with  exact  results  might  be  expected.  For  instance,  Fig  23  shows  the  variation  with 
time  of  several  response  quantities  for  system  II  of  Table  1,  when  the  normalized  excitation  frequency  is  0.5. 
Whereas  x,  and  x2  in  Figs,  7  and  8  are  essentially  sinusoidal,  x2  in  Fig.  23  is  not,  because  the  natural 


frequency  of  the  M2-K2  system  (ft2  =  1 .29)  is  close  to  the  frequency  of  the  second  Fourier  component  of  the 
friction  force  (co  =  3  x  0.5  =  1 .5),  which,  as  a  consequence,  has  a  significant  participation  in  the  motion  of  M2. 
It  is  clear  that  the  approximate  procedures  for  the  dynamic  analysis  proposed  in  this  section,  which  include 
the  first  harmonic  component  of  the  response  only,  cannot  be  expected  to  yield  accurate  results  in  cases  such 
as  this 

5  COMPARISON  OF  RESULTS  OF  DIFFERENT  METHODS 

In  order  to  assess  the  accuracy  and  applicability  of  the  different  methods  presented  in  this  paper,  results  of 
the  approximate  methods  are  compared  in  this  chapter  with  tne  corresponding  exact  results.  Values  of  the 
parameters  defining  the  system  and  the  excitation  forces  selected  to  perform  the  comparisons  are  presented 
in  Table  1.  Cases  of  pure  slip,  partial  slip  and  stuck  conditions  are  included  in  this  Table.  The  results  are 
graphically  presented  in  Figs.  24  through  38. 

The  exact  displacements  of  the  two  masses  obtained  for  Case  III  are  presented  in  Fig.  24,  in  which  they 
are  identified  by  special  symbols  defined  in  the  figure.  Results  obtained  for  this  case  with  the  harmonic 
oa.ance  method  are  indicated  by  a  dashed  line  and  those  corresponding  to  the  parametric  modification  of  the 
harmonic  balance  method  arj  shown  by  a  continuous  line.  Good  agreement  between  the  exact  values  and 
those  obtained  with  the  parametric  method  can  be  appreciated  for  high  values  of  the  parameter  p.N/Q,  (5  to 
25)  In  these  cases  the  harmonic  balance  method  is  less  accurate.  On  the  other  hand,  for  low  values  of 
pN'Q,  (1.7  and  2.5)  the  harmonic  balance  method  yields  more  accurate  results  than  its  parametric 
modification  However  the  errors  of  this  method  are  moderate  in  any  case. 

For  Cases  IV  to  XVII  of  Table  i  we  have  obtained  the  exact  values  of  the  displacements  of  both  masses 
and  the  corresponding  approximate  values  with  the  parametric  method.  Approximate  results  from  other 
proposed  modifications  of  the  harmonic  balance  method  are  in  all  cases  very  similar  to  those  presented  here. 
By  examining  these  figures,  it  can  be  concluded  that  the  parametric  method  provides  good  approximations  to 
the  exact  results,  with  the  exception  of  some  values  in  the  vicinity  of  the  resonant  frequency  of  the  M2-K2 
system  (Q2=  0.577)  in  cases  X  and  XI  (Figs.  31  and  32).  The  reason  is  that  the  natural  frequency  of  the 
M,-K,  system  is  in  these  cases  exactly  three  times  the  previous  frequency  (ft,  =  3  x  0.577  =  1.732),  and  as  a 
consequence  the  third  Fourier  component  of  the  friction  force  has  a  significant  participation  in  the  response. 
A  similar  situation  can  be  appreciated  for  Case  II  in  time  histories  depicted  in  Fig.  23,  which  were  obtained 
with  numerical  integration.  It  is  evident  that  the  higher  frequency  dominates  in  x2. 

In  order  to  assess  the  range  of  applicability  of  the  approximate  methods  we  have  examined  closely  two 
particular  cases.  Results  from  all  the  methods  considered  in  this  report  for  case  I,  with  co  =  1,  are  presented 
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m  Table  2  It  can  be  appreciated  that  the  modified  parametric  methoc  provides  approximate  values  vei  / 
ciose  to  the  exact  ones  in  all  cases.  The  c,  ladratic  method  gives  almost  the  same  results  as  the  parametric 
method.  The  errors  of  the  linear  modification  are  slightly  greater  than  those  of  the  previous  methods  The 
larger  discrepancies  with  exact  results  correspond  to  the  classic  harmonic  balance  method,  the  worst  case 
being  tiiat  in  which  pN/Q,  is  close  to  jt/4  times  the  maximum  value  of  that  parameter  (10/12.96  =  0.772).  We 
.an  ooseive  that  the  accuracy  of  the  approximate  results  in  Table  2  is  in  agreement  with  the  precision  found 
in  estimating  the  first  Fourier  coefficient  of  the  tangential  force  in  the  joint  resulting  from  all  the  methods  for 
the  same  values  of  the  system  and  excitation  parameters,  as  shown  in  Fig.  20. 

The  same  type  of  results  are  presented  in  Table  3  and  Fig.  19  for  ro  =  0.7071.  In  this  case,  the  harmonic 
balance  method  yields  the  most  accurate  results  for  all  the  values  of  pN/Qr  Indeed,  the  modified  methods 
provide  acceptable  results  only  for  low  values  of  pN/Qr  By  examining  Fig.  19,  we  conclude  that  the  reason  is 
that  in  this  particular  case  the  joint  is  actually  always  sliding  for  almost  all  the  values  of  pN/Qv  conforming  the 
assumptions  of  the  classical  harmonic  balance  method  (see  section  4.1). 

in  general,  the  modified  methods  yield  accurate  results  in  cases  in  which  the  slip-to-stuck  (r, )  and  the 
some-slip  to  pure-slip  (r2)  limits  of  pN/Q,  (see  section  3.3  and  Fig.  3)  are  far  from  each  other.  In  the  opposite 
case,  the  original  method  will  provide  better  results.  This  emphasizes  the  importance  of  calculating  those 
limits,  as  was  done  to  construct  Fig.  3  From  this  figure  it  is  possible  to  affirm  that  the  modified  methods  will 
be  more  accurate  in  the  vicinity  of  the  resonant  frequency  of  the  complete  system  (co  =  1)  where  r2/r1  =  0.07 
The  classic  method  will  be  better  when  the  mentioned  limits  have  close  values,  for  instance,  for  w  =  0.7071, 
where  ryr1  =  0.76. 

6  CONCLUDING  REMARKS 

Several  methodologies  to  study  the  steady-state  vibrations  of  the  two-body  system  depicted  in  Fig.  1  have 
been  presented.  First,  an  efficient  analytical  exact  solution  was  derived,  which  allows  to  calculate  the  actual 
values  and  characteristics  of  any  response  of  interest.  This  method  was  used  to  obtain  detailed  results  for  a 
representative  system. 

The  exact  results  show  that  the  ratio  of  the  normal  force  in  the  jo. fit  to  the  amplitude  of  the  excitation  forces, 
measured  by  the  dimensionless  parameter  pN/Q,  ,  is  a  critical  factor  in  the  system  response.  For  high  values 
of  this  ratio  the  joint  will  be  locked  in  most  cases  and  the  maximum  values  of  the  response  will  occur  for  the 
resonant  frequency  of  the  system  vibrating  as  a  unit  («=1  in  the  normalized  system).  However,  if  uN/Q,  is 
small  enough  to  allow  some  slipping  in  the  joint,  the  corresponding  energy  dissipated  will  produce  a 
significant  reduction  of  the  response  amplitude. 


For  low  values  of  the  ratio  pN/Q,,  the  joint  will  be  always  sliding,  and  the  peak  responses  will  correspond  to 
the  natural  frequencies  of  each  mass-spring  system  vibrating  independently  of  the  other  one.  Again,  even 
small  values  of  forces  in  the  joint  will  produce  significant  energy  dissipation  and  a  consequent  reduction  of  the 
response. 

An  optimum  value  of  the  parameter  pN/Q,  can  be  calculated  such  that  the  maximum  value  of  the  response 
**  interest  is  minimum,  considering  all  the  relevant  frequencies.  In  general,  the  optimum  pN/Q,  will  be  small 
enough  to  reduce  significantly  the  response  at  the  resonant  frequency  of  the  whole  system,  and  high  enough 
o  yield  appreciable  reduction  of  the  response  at  the  resonat  frequencies  of  the  bodies  vibrating 
independently. 

The  analyli'il  exact  solution  provides  the  means  to  compute  the  values  of  the  ratios  pN/Q,  corresponding 
to  the  stick-slip  transition  in  the  joint.  If  pN/Q,  is  higher  or  equal  to  this  limit,  the  joint  will  be  locked,  and  slip 
can  occur  only  for  loads  lower  than  the  limit.  Also  the  value  of  pN/Q1  separating  the  region  of  pure-slip  and 
partia1  slip  can  be  calculated.  The  compita^on  of  these  limits  is  important  to  obtain  mathematical  solutions 
consistent  with  the  physical  phenomena,  and  tc  C  'ermine  the  applicability  of  approximate  procedures. 

A  comp'-'x  variable  version  of  the  harmonic  balance  method  also  has  been  presented.  It  constitutes  a 
s  V-  wa;  *or  deriving  the  algebraic  real  equations  governing  the  problem,  and  seems  to  be  the  most 
appropiate  pproa.h  to  address  multi-degree-of-freedom  or  multi-joint  problems.  An  important  feature  of  our 
formulation  is  that  a  criterion  is  proposed  to  exactly  calculate  the  normal  load  corresponding  to  the  slip-stick 
limit  in  the  joint.  Physically  acceptable  solutions  are  possible  only  for  normal  loads  lower  than  this  limit. 

The  comparisons  of  exact  and  approximate  results  show  that  methods  based  on  considering  only  the  first 
Fourier  component  of  the  tangential  force  in  the  friction  joint,  are  very  accurate,  provided  thtw  ihe  participation 
of  higher  Fourier  components  in  the  response  is  negligible.  Such  a  participation  can  be  prevented  by 
designing  systems  in  which  the  higher  frequencies  are  not  odd  multiples  of  the  lower  ones. 

Modifications  to  the  harmonic  balance  method  are  also  proposed.  They  consist  in  estimating  the  first 
Fourier  coefficient  of  the  tangential  force  in  the  joint  with  approximate  heuristic  rules.  The  improvement 
achieved  in  cases  of  partial  slip  with  these  modifications  shows  the  relevance  of  a  proper  estimation  of  the 
first  Fourier  coefficient  of  the  friction  force.  In  cases  in  which  the  joint  is  sliding  during  an  entire  cycle  the 
classical  assumption  that  such  a  coefficient  is  equal  to  4/ji  will  provide  the  best  results.  If  the  slipping  is  only 
partial,  the  rules  proposed  in  the  modified  harmonic  balance  methods  will  produce  better  results.  The 
quadratic  rule  seems  to  be,  in  these  cases,  a  good  compromise  between  simplicity  and  accuracy. 

The  agreement  of  most  of  the  exact  results  with  the  approximate  ones  supports  the  assumption  that  the 
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actual  response  of  the  two-body  system  studied  in  this  work  is  nearly  harmonic.  This  behavior  can  also  be 
expected  in  multi-joint  and  multi-body  systems,  provided  the  effect  of  high  frequencies  is  reduced  to  a 
minimum  by  means  of  an  appropiate  design  of  the  system.  This  would  allow  the  application  of  the 
approximate  methods  presented  here  to  the  analysis  of  such  systems,  resulting  in  a  substantial  simplification 
of  the  mathematical  problem. 
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COMPARISON  OF  RESULTS  FROM  DIFFERENT  METHODS 


xt  =  normalized  displacement  of  mass  M1  of  Fig  1. 
x2  =  normalized  displacement  of  mass  M2  of  Fig.  2. 
E  =  energy  dissipated  per  cycle 


C,  =  First  Fourier  coefficient  of  the  friction  force 
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Fig.  6:  Maximum  Exact  Displacement  of  Mass  M2  of  Fig.  1. 
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Fig.  7:  Exact  Rasponsa  Time  Histories  for  a  Casa  with  tha  Joint  Stuck 
(Casa  I  of  Tabia  1.  ,^(1^12.96.  a. =1.0) 
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Fig.  8:  Exact  Response  Time  Histories  for  a  Partial-Slip  Case 
(Case  I  of  Table  1.  ^N/Q  =6.91.  *.  =  1.0) 
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Exact  Response  Time  Histories  for  a  Pure-Slip  Case. 
(Case  I  of  Table  1,  //N/Q^O.67,  u> =1.41) 
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Fig-  10:  Discontinuities  in  Response  Functions  in  a  Case  of 
Partial  Slip  in  the  Joint. 
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Fig.  13:  Energy  Dissipated  per  Cycle.  Csse  i  of  Table  1.  <#>0.7071. 
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Fig.  20:  First  Fourier  Coefficient  of  the  Friction  Force. 
Case  I  of  Table  1,  u>  =  1.0. 
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Fig.  21:  Variation  of  Cf  with  D  in  the  Approximate  Harmonic 
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Fig.  27:  Comparison  of  Exact  and  Approximate  Results.  Case  VI  of 

Table  1. 
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Fig.  26:  Comparison  of  Exact  and  Approximate  Results.  Casa  VII  of 

Table  1. 
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Fig.  32:  Comparison  of  Exact  and  Approximate  Results.  Case  XI  of 

Table  1. 
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Fig.  33:  Comparison  of  Exact  and  Approximate  Results.  Case  XII  of 

Table  Y. 
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APPENDIX  A 


cos(^  +  ft)  =  -sin  ft  sin  f  +  cos  ft  cos  + 


Substitution  of  y  and  y2  into  1.6  yields 
g,(0,^)  =  [Btcos^  -  Dtsin^]o> 


g2tO.^)  =  [!B2cos/$  -  DjSin/?)  cos^  -  (B2sin/?  +  D  cos  ft)  sin^]a/ 


Define 

q,  =  [B1  -  (BjCos^  -  DjSinyS)]^ 
p,  =  [~D)  +  (B2sin/S  +  D2cos  ft)~]w 
From  1.1  to  1.4  we  can  also  define 
b„  =  f,(0)  =  O, 
b12  =  f2(0)  =  -c,o, 
b,3  =  f,(0)  =  o; 

bu  =  V°>  = 

Expression  (14)  follows  from  substitution  of  A. 7  to  A 12  into  expression 
(13). 


c)  Functions  defining  the  solution  of  equation  (8) 

The  following  functions  appear  in  expresion  (17) 

f 5(t)  =  exp(-Ct)  sin(\/T-CJ  t) 

f#(t)  =  exp(-Ct)  cos<yT-CJ  t) 

g3(t.^)  =  2C<i>[sin(wt+^)  +  sinUt+^  +  .?)]/L  +  (1-uiJ)[cosM+^)  + 

cos(<«t+^+^)]/L 


where 


C  =  C,  /(Win -m)  +  C2  s/ km 


L  =  (1  -  w2)2  +  (2Cw)2 


APPENDIX  B.  Derivation  of  equations  (22)  to  (24) 


Equation  (22)  is  obtained  by  subtracting  (20)  from  (19).  From  (21).  E2  can  be 
expressed  as 

E,=  [ A , f , (t , )  +  A2f2(t,»  +  g,(t,  f)  -  g3(t,.  f)  -  E,f ,9(t1)]/fe(t1) 

Substituting  this  result  and  equation  (22)  into  (19)  and  solving  for  Er  we 
obtain 

E  =  A,f7(t,)  +  Aifa(ti)  +  A3f9(t1)  ♦  A4f10(t,)  +  g4  (t,.  f)  +  C,fQ 
where  the  following  definitions  have  been  adopted 


W  =  (1-k)f2(t,)  -  f^t^f.tt,)/^,) 


f9(t,)  =  kf3(t,( 


fJV  =  kW 

g4(trf)  =  (1-k)g|(t|^)-g1(t1^)fett1)/f#(t|)-gJ(tJ^)+gJ(t1^)fa(t1)/fe(t1)+kg2(t1^) 
c,  =  ( 1  — k)/(  1  — m)  +  k/m 


Backsubstitution  of  Ef  into  E}  yields  equation  (24)  with  the  following 
de.  nitions 


fn<V  =  tW  -  WW»W 


=  cf2(ti'  -  WW^'W 


f,3(ti>  =  -  WW'W 


=  -  ViWW 
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B.l 


B.2 


g5(t,'  **>  =  <v  -  93(t,.  -  g4(t,.  j»)fs(t)]/f,(tt) 

C3(t,)  =  - 

Note  that  g4(t(,  and  gs(tr  f)  are  eventually  linear  combinations  of  the 

harmonic  functions  g(t.  f  )  (j=  1,  2,  3)  defined  in  Appendix  A.  Consequently 
they  can  be  expressed  in  the  form  piysin^  +  gtt^cos^  that  is  used  in 
equation  (29). 
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